Radon transform and kinetic equations in tomographic representation 
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Abstract 

Statistical properties of classical random process are considered in tomographic representation. The 
Radon integral transform is used to construct the tomographic form of kinetic equations. Relation of 
probability density on phase space for classical systems with tomographic probability distributions is 
elucidated. Examples of simple kinetic equations like Liouville equations for one and many particles 
are studied in detail. 



O 1 Introduction 



The kinetic equations describing the behavior of classical system can be obtained using procedure sug- 
gested by Bogolyubov [I] where the chain of connected equations corresponding to N and N +1 particles 
in the system was studied (see also |2J). The simplest Liouville equation for N + 1 particles probabil- 
• i-h [ ity density on phase-space provided the more complicated equation for N particles probability density 
^ ' after averaging over (N + l)th particle degree of freedom. The appropriate reduction procedure of the 
■ Bogolyubov chain of equations gives other kinetic equations, for example Boltzman equation with the 
collision term. The probability density on phase-space f(q,p, t) = f(q\,q2,---,qN,Pi,P2,---,PN,t) depends 
on positions and momenta of the particles. Recently it was suggested to use for description of quantum [3j 
and classical [I] states the tomographic probability distribution function (state tomogram). For classical 
particles the tomogram is Radon transform [5] of the probability density f(q,p, t) on the phase-space. The 
advantage of the tomographic probability representation is connected with the fact that in this represen- 
tation both classical and quantum states are described by the identical objects-tomographic probability 
distributions. The physical meaning of the tomographic probability distribution is the following one. For 
one particle it is distribution of the particle position X but this position in measured is the reference 
frame in the system phase-space which rotated by angel 9 and before the rotation the axes q — » sq, 
p — ► s _1 p where s it the scaling parameter. Then the scaled and rotated frame can be parameterized 
by two real parameters /x = scos# and v = s _1 sin^. For one particle the tomogram is the probability 
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distribution function w{X,p,v) of random position X and two real reference-frame parameters fi, v. 
For N particles the tomogram is the function w(X,j2,P) where X is vector with particles positions Xj„ 
j = 1,2, ...,N and fl and v vector components [ij = sjcosOj and Vj = s~ 1 sin0j describe the scaling 
and rotation parameters of reference-frame axes associated with jth degree of freedom. The aim of this 
work is the study of classical kinetic equations in tomographic probability representation. We focuse 
on the Bogolyubov chain written for the function f(q,p,t) and apply the Radon transform to write the 
chain of kinetic equations. On the other hand the consideration can be extended to quantum kinetic 
equations since in the tomographic picture the connection of tomograms for N + 1 and N particles is 
identical for classical and for quantum domains. This fact corresponds to analogous relations for Wigner 
functions [6] in quantum domain but the Wigner function can take negative values and it means that the 
function is not probability distribution function. Some aspects of the tomographic description of classical 
and quantum states were considered in [7] , [8] , [9] , [10] , [11] , [12] , [13] , [2] . The paper is organized as follows. 
In Sec. 2 we review the approach with the tomographic Radon map of classical probability distribution 
both for one particle and for N particles. In Sec. 3 we obtain the Liouville equation in the tomographic 
probability representation. In Sec. 4 we review the procedure with averaging the Liouville equation used 
in Bogolyubov chain approach. Also we consider the Radon representation. In Sec. 5 the perspectives 
and conclusions are given. 



2 Probability distribution and tomograms 

Let us consider a particle with fluctuating momentum p and position q. The state of such particle is 
associated with probability distributions f(q,p) on the particle phase space. The probability distribution 
is nonnegative and normalized function, i.e. 

f(q,p) > o (i) 

and 

J f(q,p)dqdp=l. (2) 

The function f(q,p) can be considered also as generalized function, e.g. Dirac delta-function. For many 
particles the state of the system with fluctuations is described by the probability distributions f(q,p) 
where q = (q±, (fe, • ■> qN), P = (pi,P2, --,Pn) but this function is nonnegative and normalized as in one- 
dimensional case. The Radon transform of the probability distribution function f(q,p) is gives by the 
formula 

w(X, = J f(q, p)5{X - /j,q - vp)dqdp (3) 
where the arguments X,/i,i/ are real numbers. The function w{X 1 ^ 1 v) is nonnegative and normalized, 



i.e. 



Jw(X,fi,u)dX = 1. (4) 

This can be checked using integrations over X in Q and property of Dirac delta-function. The function 
w(X, fx, u) is called symplectic tomogram. The physical meaning of the symplectic tomogram is the 
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following one. If the parameters n and v are represented in the form [i = scosO, v = s -1 sin# the 
symplectic tomogram is the probability distribution of the particle position measured in the reference 
frame in the particle phase-space with following properties. Let us first scale the axes of the reference 
frame q — > sq = q', p — > s~ l p = p' and than rotate new resealed axes q' and p 1 by the angle 6. Then we 
get reference frame with axes where q" = q 1 cos9 + p'sin9 and q" = —q'sm.9 + p' cos6. The coordinate 
X it the particle position measured in the reference frame with new coordinates q" and p" . Thus taking 
H = s cos 6, v = s -1 sin 6 we have 

X = /j,q + up. (5) 

The Radon transform ([3]) has inverse, i.e. 

f(q,p) = A / M, v)e l{X - m - vv) dXdp,dv. (6) 



4vr 2 

The symplectic tomogram has the homogeneity property 

w(XX, X/i, \v) = r^MX, /I, v). (7) 

This property follows from homogeneity property of Dirac delta-function 5(\y) = ^5(y). For many 
particles one has the multidimensional Radon transform 

f n 

w(X, fl,v)= I f{q,p) 5(X k - pL k q k - v k p k )dqdp. (8) 
J k=X 

This symplectic tomogram also determines the probability distribution f(q,p) on the phase-space of the 
system, i.e. 

1 r N 

f(q,p) = j^yj J w(X,M(Yl e^-^'^)dXdfldu. (9) 

The symplectic tomogram ([8]) is nonnegative and normalized, i.e. 

w(X,fl,u) > (10) 

and 

' w(X,ft,P)dX = 1. (11) 



The physical meaning of the symplectic tomogram w{X,jl ) V) is the following one. It is the joint proba- 
bility distribution function of N random variables (positions) X k measured in specific reference frame in 
the system phase-space subjected to operation of rescaling q k — > s k qk, Pk —> p k and then rotated by 
the angle Q k . Then the position X k is measured in the reference frame where all the axes are transformed 
using this prescription. Important property which follows from the physical meaning of the tomogram 
w(X, p, v) is the reduction property, i.e. 

w(X, p, v)dX N = w(X', p', if). (12) 
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Here X' = {X\, X 2 , X^-i),^ = (/xi, /i2, P-n-i) and z/ = z/ 2 , i/jv-i). This reduction property 
corresponds to the reduction property of the probability distribution function 

y f(q,P)dq N dp N = f((f,j?) (13) 

where <f = (g, q^-i), p* = (p, ...,Pn-i)- The function f(q^,p f ) is the probability distribution for 
the subsystem of N — 1 particles. This probability distribution determines the symplectic tomogram 
w(X', //, z7) by using the Radon transform 

- N-l 

/ II S ( X * - VkQk - VkPk)dqdp = w(X', jt, if). (14) 

J k=i 

The symplectic tomogram of multidimensional system has the homogeneity property 

N 1 

w(\iXi, \ 2 X 2 , \ N X N , \ 2 fi 2 , \ N fl N , \iV 1: \ 2 U 2 , \ N V N ) = (]J JVl") W ( X > & (15) 

fc=l l Afe l 

The function must satisfy the condition 

at 

/ /x, i?) J] e l{Xk -^ qk - UkPk) )dXdpdv > (16) 
fc=i 

since this integral determines the probability distribution of the system on the phase-space. 

3 Liouville equation in tomographic representation 

The simplest kinetic equation in classical statistical mechanics is Liouville equation for the probability 
distribution f(q,p,t) of one particle on the phase-space. This equation reads 

df(q,P,t) df(q, p, t) dU (g) dfjq, p,t) _ Q 
dt dq dq dp 

Here U (q) is potential energy and this equation corresponds to the Newton equation of motion 
In the case of many degrees of freedom the Hamiltonian reads 

N p 2 

H(q l ,q 2 ,...,q N ,Pi,P2, -,Pn) = ^2 ~^ + U(qi,q 2 , ...,q N )- (19) 
The Newton equation of motion corresponding to this Hamiltonian has the form 

A = (20) 
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The Liouville kinetic equation in this case reads 

d/(gi,g2, -,QN,Pi,P2, -,PN,t) A df(qi,q 2 , ...,q N ,pi,p 2 , -,PN,t) 

at + ^ Pi o qj 

df(q 1 ,q 2 ,...,q N ,p 1 ,p 2 ,...,p N ,t)dU(q 1 ,q 2 ,...,q N ,p 1 ,p 2 ,...,p N ) = Q ^ 
p{ dPi <% 

The Liouville kinetic equation can be rewritten in tomographic probability representation. It means that 
we replace the probability density f(qi, q 2 , qN,Pi,P2, ■■■,PN,t) by its Radon transform and write down 
equation for the tomographic probability density. For one degree of freedom we use the correspondence 
rules 

f(q,p,t) < — > w{X,n,v,t); 

d d 

pf(q,p,t) < — > -(—y 1 —w(X,fi,v,t); 

d d 
■7^f{q,P,t) < — > n—w(X,fi,u,t); 

d d 

qf(q,p,t) < — ► -(—)' 1 —w{x,fi,u,t); 

^/(g,P, t) <— > v^w{X, fi, v, t). (22) 

Then we get the Liouville kinetic equation in the tomographic probability representation 

dw(X, fj,,u,t) dw(X, fj,,i/,i) 
dt ^ dv 

dq [q [ dx> V dx " u - [26) 

For N degrees of freedom the correspondence rules have the form 

/(gl,g2, ...,qN,Pl,P2, ...,PN,t) < > w(Xi,X 2 , ...,X N ,fJ,i,fl 2 , -,fJiN, "I, "2, V N ,t)] 

d d 

Pjf(qi,q2, ...,qN,Pi,P2, -,PN,t) < — > ~(^") _] l ~frJ~ w ( Xl - • X2 ' -, X N, 1^1,1^2, -,m, v l,V2, ...,u N ,t); 

d d 
^-/(gi,g2,-,giV,Pi,P2, -,PN,t) < — > n j -^-w(X 1 ,X 2 , ...,X N ,n 1 ,H2,...,n N ,v 1 ,v 2 ,...,v N ,t); 

d d 

gj/(gl,g2, -,qN,Pl,P2, -,PN,t) < > -(——y 1 —w(Xi,X 2 , ...,X N , 1^1,1*2, -, f*N, "l,"2, ...,l>N,t); 

OJLj Oflj 

d d 

7^-/(gl,g2, ...,qN,Pl,P2, -,PN,t) < > V j Q^-w(X 1 ,X 2 , ...,X N , m, /J.2, UN, Vl,V2, -,VN,t). (24) 

The operator (^) _1 acting on a function f(X) gives the function $(X) satisfying the equality 
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In case of representing the functions ip(X) and &(X) in form of Fourier integrals, i.e. 

ip(X) = J y>(k)e ikX dk (26) 

and 

fcpf) = J 4>(k)e ikX dk (27) 
we define the action of the operator (^) _1 on the function tp(X) by equality 

( ^ rV(x) = / y>^y kXdk ( 2§ ) 

which removes the ambiguity in choice of constant in the function &(X). Now we can write down the 
system with many degrees of freedom in the tomographic probability representation. To do this we make 
Radon transform in (I2ip . Then we get 

dw{X,p,,v,t) ^ dw(X,jl,u,t) 

E v dU , d i d . d 1 d 



' x dqj dXi dm ' dX 2 dfi 2 



, d ! d dw(X,n,v,t) 

oXn OfXN 8Xj 

The written equation describes evolution of tomographic probability distribution function which corre- 
sponds to standard Liouville equation for the probability density on the phase-space. 



4 Reduced Liouville equation 

The Bogolyubov chain of equations for probability distribution function of one particle can be obtained 
using Liouville equation for many particles. The ansatz is the following. One integrates the Liouville 
equation over N — l pairs of variables q 2 , (ft, qN,P2,P3, ■■■■>PN- Then one has the equation for one degree 
of freedom only with the contribution of the interaction of one particle with the rest. We provide this 
derivation considering system of two particles with the Liouville equation of the form 



d ,, , df(q 1 ,q 2 ,pi,p 2 ,t) df(q 1 ,q 2 ,pi,p 2 ,t) 

■Klf{QUQ2,Pl,P2,t) +P1 ^ hj>2 ~ 

ot oqi oq 2 

dU{qi,q2) df(qi,q 2 ,pi,P2,t) _ dU(qi,q 2 ) df(qi,q 2 ,p 1 ,p 2 ,t) 
dqi dpi dq 2 dp 2 



0. (30) 



We integrate this equation over variables of second particle q 2 ,p 2 - Let us assume that potential energy 
depends on distance between the particles, i.e. 

U( qi ,q 2 ) = U(\ qi -q 2 \) = U(\X lt2 \), 

Xi, 2 = q\ - q 2 . (31) 



(3 



Then we get 



a ~ a ~ r a 

-^f{qi,pi,t) +Pig^f(qi,Pi,t) + J P2-g^f(qi,q2,Pi,P2,t)dq 2 dp 2 - 

f , d d 

-J U\\X h2 \){—\qi - q 2 \)-^f(qi,q 2 ,Pi,p 2 ,t)dq 2 dp 2 + 

+ f U'(\X lt2 \)(^-\q 1 -q 2 \)^-f{q 1 ,q 2 ,p 1 ,p 2 ,t)dq 2 dp 2 = 0. (32) 
J oqi dp 2 



f(qi,Pi,t) = J f(qi,q2,Pi,P2,t)dq 2 dp 2 . (33) 



We use the notation 



This equation can be transformed and analogous procedure can be done with Liouville equations (I21j) . 
Now we make the reduction of tomographic Liouville equation (|29|) on the example of two particles. The 
equation (j29j) for two particles reads 

dw(X 1 ,X 2 ,^,fi 2 ,u 1 ,iy 2 ,t) , d d 

di du^ + n 2 -^)w{X 1 ,X 2 ,fi 1 ,fi 2 ,iy 1 ,v 2 ,t) - 

(34) 

We can get the reduction of this equation for tomogram describing the state of one particle integrating 
it over variable X 2 . One has 

dw(Xi,m,vi : t) d 

-Q t /ii— u>(Xi,/ii, ,i/i,t) + 



The obtained equation (|35|) is consistent with taking Radon transform of the reduced Liouville equation 
(|32|) . The Bogolyubov chain of equations in case of many particles can be obtained by analogous procedure 
applied to equation ([29]) . In this case we integrate the equation over variables X 2 ,X%, ...,X^f. 



5 Conclusion 

We point out the main results of the work. We considered Liouville equation for N particles and using 
the integral Radon transform of the probability density on phase-space rewrote the equation in the 
tomographic probability representation. The Bogolyubov chain of the equations obtained by averaging 
the Liouville equations over subsystem degrees of freedom is also written in the tomographic probability 
representation. We hope to extend in future publications the tomographic representation analysis to 
the quantum kinetic equations using the Bogolyubov chain approach to Moyal equation [15] for Wigner 
function. 
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